Full counting statistics of spin transfer through the Kondo dot 
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We calculate the spin current distribution function for a Kondo dot in two different regimes. 
In the exactly solvable Toulouse limit the linear response, zero temperature statistics of the spin 
transfer is trinomial, such that all the odd moments vanish and the even moments follow a binomial 
distribution. On the contrary, the corresponding spin-resolved distribution turns out to be binomial. 
The combined spin and charge statistics is also determined. In particular, we find that in the case 
of a finite magnetic field or an asymmetric junction the spin and charge measurements become 
statistically dependent. Furthermore, we analyzed the spin counting statistics of a generic Kondo 
dot at and around the strong-coupling fixed point (the unitary limit). Comparing these results with 
the Toulouse limit calculation we determine which features of the latter are generic and which ones 
are artifacts of the spin symmetry breaking. 

PACS numbers: 72.10.Fk, 72.25.Mk, 73.63.-b 



I. INTRODUCTION 

It has been known for some time that the current 
noise spectra of mesoscopic systems can reveal informa- 
tion about transport processes which is not accessible 
via measurement of the average current alone. Indeed, 
it has been shown by Schottkj^ that the characteristics 
of charge carriers can be investigated by measuring noise 
spectra. The noise is the second moment of the distri- 
bution function which gives the probability that a cer- 
tain amount of charge is transported within a certain 
time interval T. Recent developments in experimental 
mesoscopic physics have facilitated the measurement of 
even higher order moments of this distribution,^ There- 
fore, it is desirable to calculate the complete statistics, 
which is often referred to as the full counting statistics 
(FCS) . The FCS has been found for the charge transport 
in non-interacting systems by Levitov and Lesovife^- and, 
currently, interaction effects are under intensive investi- 
gation. 

Recent years have witnessed a soaring interest in the 
field of spintronicSf^ where the principal objective is to 
gain a level of control over spin transport (which com- 
prises creation of current and its measurement) compa- 
rable to that over charge transport. Indeed, the usage 
of the electron spin degree of freedom offers many entic- 
ing opportunities which are currently being investigated. 
Therefore, in order to gain a thorough understanding 
of spin-related processes taking place in mesoscopic sys- 
tems, one has to understand the properties of spin current 
to the same level of accuracy as the charge current. 

A method of measuring spin currents and some intrigu- 
ing consequences of consecutive measurements of differ- 
ent spin components have been presented by Di Lorenzo 
and Nazarov."'' Although some general statements about 
spin transfer FCS have been made there, to the best of 
our knowledge, there has so far been no derivation of 



the spin transfer FCS for a concrete system. The aim 
of this article is to fill this gap. While the charge FCS 
was shown to be binomial, at least in the universal linear 
response, zero temperature limit f^ it is not immediately 
clear what the nature of the spin statistics should be: the 
natural guess that it is the inverse binomial distribution^ 
(due to bosonic nature of spin excitations) turns out to 
be wrong. The purpose of this article is to investigate 
the spin transfer FCS for the transport through a Kondo 
dot, where non-zero spin currents can be relatively eas- 
ily generated by applying a finite magnetic field and a 
bias voltage. Furthermore, we go one step further and 
calculate the combined FCS of charge and spin transfer, 
which will allow us to examine cross-correlations of these 
two currents and to obtain spin-resolved charge transfer 
statistics. 

We shall perform this analysis on two distinct systems 
which represent two aspects of the Kondo physics. First, 
we will investigate a Kondo dot at the Toulouse point, 
where for a certain choice of coupling constants, the sys- 
tem can be diagonalized exactly by means of bosoniza- 
tion and subsequent refermionization. Second, we will 
go to the strong-coupling regime where the dot spin is 
hybridized with the leads spins to form a singlet state. 
Virtual excitations of this singlet state can be treated 
perturbatively and will be shown to produce a finite spin 
current noise. This article will close with a comparison 
of these two regimes. 



II. HAMILTONIAN 

The system under investigation consists of two non- 
interacting leads and an interacting dot level. The 
left/right leads are being held at chemical potentials 
Mi/fli respectively, where V = (/i^ — ^iRj/e is the ap- 
plied voltage. The position of the bare dot level and the 



electron-electron repulsion are such that the dot is always 
occupied by exactly one electron. The result is the two- 
channel Kondo Hamiltonian,— H = Hq + Hk + Hb + ^, 
where (in units oi e = h = vp = 1), 

a=LRa=n 
aP=LR v=x,y,z 

-^iBghT"" =: -At' , 
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In this equation, r^'^'^ are proportional to the Pauli ma- 
trices describing the impurity spin, J"^ for a, f3 — L, R 
and v = x,y, z are coupling constants for the possible 
transport processes, h is the magnetic field applied to 
the dot level. Furthermore, 
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(2) 



for a, f3 = L,R are the components of the electron spin 
density at the impurity site. 

The aim of FCS is the calculation of the cumulant gen- 
erating function (CGF) x(/i. A) for the probability dis- 
tribution function P{s, q) to transfer q units of charge 
and s units of spin during the waiting time T. In our 
case, as we are interested in the combined statistics of 
spin transfer and charge transfer, we define x(Mi ^) — 
X^s g e*'*''e**'^P(s, q). Partial derivatives of this function 
then allow us to calculate arbitrary correlation functions 
of charge and spin transport and cross-correlations be- 
tween both. 

In order to calculate spin and charge current in this 
system, we proceed as usually in tunneling setups and 
define charge and spin flavor number operators 

^c = 2 E E " / dx^\Ax)i'cca{x) , (3) 

The time derivatives of these operators are proportional 
to the charge and spin currents, respectively. Hence, the 
corresponding charge and spin currents can be calculated 
by means of the Heisenberg equation Ic^s — i[H,Nc,s]- 
The general approach to calculate higher order cumu- 
lants of these operators involves including fictitious mea- 
surement devices for both charge and spin current in 
the Hamiltonian.^ This leads to additional terms H„i — 
fils + XI c- As in our case the particle number opera- 
tors Nc and Ns commute, both terms can be removed 
by a unitary transformation. The counting fields /^ and 
A then appear as exponentials in those components of 
the Hamiltonian which transfer spin or charge in either 



directioujiS from the left lead to the right, T^'*, and back- 
wards, T^'*'. As in this case the underlying processes for 
spin and charge transport will be different, the tunnel- 
ing Hamiltonian separates into spin and charge parts and 
becomes [see Eq. (ITO)) ] 
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(5) 



where piit) and A(i) are explicitly time-dependent func- 
tions defined on the Keldysh contour C . The count- 
ing fields are only "switched on" during the measure- 
ment time T, such that we can define them as /i(t±) = 
±fie{t)e{T - t) and analogously for \{t±). The CGF is 
then given by the expectation value^- 



x(m, A) = (Tc 



-I J^ dsT''^ (s) 



(6) 



where Tc denotes the contour ordering operator. In or- 
der to calculate this function, we employ the general- 
ized Green's function formalism that has previously been 
used successfully to calculate the FCS of other interact- 
ing systemsi-Siii We assume A(s) and ^(s) to be arbitrary 
functions defined on C that have slowly varying Keldysh 
components \±{t) and /i±(t), respectively. Assuming T 
to be large, we can neglect processes which are due to the 
switching of the counting fields, and we can introduce the 
adiabatic potential U{fi±, X±) which entails 



x[^^±it),x±it)] ^e-*^o'^''*^[^±(*)'^±wi 



(7) 



from which the statistics can be recovered using the equa- 
tion lnx(/i. A) = —iTU{fi, A) as C/(/i, A) will turn out to 
be time- independent. By comparison of ([6]) and ([7]), one 
notices that the adiabatic potential can be calculated due 
to the following relation (which is a consequence of the 
Feynman-Hellmann theorem) 



dU 



^a7 



(8) 



/^A 



or, equivalently, via the partial derivative with respect to 
/i_ . The expectation value with subscript /^A is defined as 
expectation value in the interaction picture with respect 
to the tunneling Hamiltonian T^"^, ie.-i^ 

{A{t))^, = -^-L— (tc {A(t)e-/c---(^)}) . (9) 

The emerging expression is slightly more complicated 
than the usual Hamiltonian formalism. However, the 
great advantage of using the derivatives of U is the emer- 
gence of Green's functions (the time dependence of which 
is dominated by T^^^) for which a Dyson equation can be 
obtained. Then, the adiabatic potential U{^, A) can be 
recovered by integrating ([5]) (or, equivalently, the deriva- 
tive with respect to /i^, the result is the same) which, in 
turn, leads immediately to the CGF x(/i. A). 



III. DYSON EQUATION AT THE 
TOULOUSE POINT 



For a general constellation of system parameters, 
x(/i, A) cannot be calculated exactly. The perturbation 
theory in the exchange couplings is well known to di- 
verge logarithmically and so is not very useful in calcu- 
lating the FCS. However, it has been realized by Toulouse 
that the Kondo Hamiltonian can be reduced to an eas- 
ily integrable shape for a special parameter setJ^ In 
our notation this so-called Toulouse point is attained for 
Jf = Jf =: jf^jLR ^ jRL and J^^ = Jf^ = 04 
The only allowed processes are then: (i) spin-flip tun- 
neling which transports one unit of charge but does not 
contribute to the spin current, and (ii) spin-exchange be- 
tween one lead and the dot which only transports one 
unit of spin with no associated charge transport. Fur- 
nishing both processes with the respective counting fields 
^,{t) and \{t) the tunneling operator becomes 
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i^l.i^Rie-'^'^T- +i,l,i;R^e-'^' 



+ i^i.i^Lie'^'^T- +i:l,4^L^e'^'\+ 



2^+ 



(10) 



In the next step, we proceed along the lines of Schiller 
and HershfieldfS First, we bosonize the Hamiltonian by 
introducing four bosonic fields describing total charge (c), 
charge flavor (f), total spin (s) and spin flavor (sf). Then, 
we perform an Emcry-Kivclson^'^ rotation and refcrmion- 
ize the Hamiltonian. After going to the Toulouse point 
J^^ ~ J^^ =: Jz '■— 2n, the total spin field decouples 
and the resulting Hamiltonian is quadratic in the three 
remaining fermionic fields. The unperturbed part reads 



Ho = 
Hb = 



» y^ / dxil)l{x)dxipi,{x) ; 

A{dU-l/2), 

V fdxij}{x)^f{x), 



(11) 



where d — ir^ is the fermion operator describing the dot 
spin and v G {c, /, sf}. The tunneling part becomes 
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'/ - e'^/^dV-/] , (12) 



where oq is the lattice constant of the underlying lattice 
model which emerges in the bosonization process. Note 
that T^'^ only contains ipf and ipsf which means that 
'0c is completely decoupled and can safely be discarded. 
Physically, this is due to the fact that the total charge 
in conserved in the system. Next, we split the relevant 
fields into hermitian components (Majorana fermions) 
according to Vs/ = (C + «'7)/\/2, 0/ = (C + ie)/V^, 
d = {a + ib)/\/2, and obtain the Hamiltonian 

H = Ho + iV I dx (:{x)e{x) - i/\ab + T^^ , (13) 



where the tunneling operator is given by 

r^^ = i J- sin(/i/2)a^ + i J+ cos(^/2)6C 
+ i J_ cos{fi/2)ari — ij^ sin(/i/2)677 
+ i J_L sin(A/2)aC + i J_L cos(A/2)ae . (14) 



In this equation, wc have introduced J± = {J^^ ± 
Jf^)/(2v%i^) and J_l = Jf ^/(\/2^). According to 
(IH), the next step is to calculate the expectation value 
of the partial derivatives of T^^^ which can easily be cal- 
culated from mi). If we calculate the derivative with 
respect to, say, A_, 

^X-/ ^ = -^[cos(A_/2)G-^- -sin(A_/2)G-r] , 

(15) 
we obtain an expression which contains the Green's func- 
tions 

G„y(i,t') = -z (ra(t)C(t')>pA (16) 

G--{t,t') = -t{Ta{t)e{t'))^^ (17) 

where the time arguments t^t' are located on the "— "- 
branch of the Keldysh contour and, consequently, are 
time-ordered. Expanding P^ to first order in J±, one 
can express it in terms of the exact dot Green's function 
Daa{t,t') — —i{Tca{t)a{t')),. Assuming the counting 
fields /i±(t) and A±(i) to be time- independent on the 
respective branches of the Keldysh contour, we have 



Gac (^) = ^J± 






^'-^aa G'j^ 



S+Daa G^^ 



(0) + - 



C+D-+G 



(0) + - 



(18) 



where s± = sin(A±/2) and c± = cos(A±/2). An anal- 
ogous equation can be derived for Gae- The exact dot 
Green's function can be derived by means of a Dyson 
equation which in Fourier space reads D = D(°)-|-DSD'°^ 
The matrix of dot Green's function is defined by 



/ Daa- D-+ D-,- D-+ \ 

\ ba ha hb 



aa ab 

Dt+ D+- 



D7+ 



bb 
^bb 



(19) 



and the self-energy S is defined as a block matrix cou- 
pling the Majorana components a and b of the dot oper- 
ator, 



and Mk 



-Mab- Next, we have for the Albt element 



Maa Mab 

Mbr, Ml 



bb 



(20) 



It turns out that this matrix only depends on A := (A_ — 
A+)/2 and fj, := {fi- — /i+)/2. We therefore introduce the 
shorthand notation s^a = sin(/i, A) and c^^a = cos(/x. A). 
The non-diagonal elements are then given by 



Mab = J-J+ 



s^G^^ 





(0)+- 






(21) 
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M, 



bb 



Ji 



r<(o) — 
„ --(0)+- 






(0) — 



(22) 



These three matrices do not involve the charge counting field A which is a consequence of the fact that the A- 
dependent part of the tunneling operator (|14p only couples to the Majorana fields a. The only A-dependent contribu- 
tion in the Dyson equation arises from the Maa component which is given by 



Maa - J- 



G 










G 



(0)- 
CC 



-caG, 



-caG(?+- - .AGi?^ 



(0)- 



5aG 



G 



(0)+4 



(0)- 



(23) 



As the zero order lead and dot Green's functions are well 
known, ^^ solving the Dyson equation amounts to invert- 
ing a 4 X 4-matrix. The ensuing calculation is lengthy 
but straightforward. The result is an expression for the 
exact dot Green's function, which can be used to evalu- 
ate (fT5|). The equation can then be integrated to obtain 
the adiabatic potential C/(/i, A) which in turn determines 
the CGF. 



IV. CUMULANT GENERATING FUNCTION 

Finally, one obtains the CGF which describes the com- 
bined statistics of charge and spin transport and which 
is the main result of this article. 



'°^<''-^>^^r^'°{'+Dfc){''(='"+' 



tively, the 'transmission coefficients' are given by 



nL{l - riR) [1 + 2nF{nF - 1)] F^Pi 

[npil - ur) + ni(l - np)] F+Fi(F_F+ + A^) 

Lo'^ [npil - nji) + nL(l - iif)] F-F^ , 

[nR{l - ul) + nL{l - ur) - 1] npinp - 1)F^F^ 



/34 



V4 



= nL(l -r^fl)n-F(l - ?i_f)F+F^ , 

= 2uj [uFiuL - nn + 1) - ul] VF-F+F^A , (26) 



-i2/j. 



-t- 



J2 ["p( 
p=± 



,i2p\ 



1) + yp{e'^ 



g-'V^gVA 



gj:(p+pA) _|_ g-i(p-pA) 



+ Pp U^^^'+P^'^ + e-'2(M-pA) _ 2 



)]} 



_ 2") where all the quantities with opposite subscripts are 
given by the same expressions with interchanged R and 
L indices. From ([M]) . one can easily calculate all correla- 
tion functions of charge and spin currents. In particular, 
by setting /i = 0, one immediately recovers the known 
charge transfer statistics for this system.-^ Before return- 
ing to the complete statistics ([M|) . we shall investigate in 

(24) niore detail the spin transfer statistics. By setting A = 0, 
we obtain the CGF for the spin transfer. 



where F^ — Jf /2 and 



D{u) = [c.2-A2-F+(Fi+F_)]- 



u\V, 



,2 /T-_^ ^ p_ 



(25) 



Defining nF{uj),nii{uj) and nL{uj) to be the Fermi func- 
tions with chemical potentials 0, —V/2 and V/2, respec- 



/•OC 7 /' 

lnx(A*) = T j —\n\l + T^a{nF-n){e^^^-e-'^') 
+ Tis [n{l - np) + nF{l - n)] {e'^' + e"*^ - 2) 
+ r2nf(l-nj^)(e2*^ + e-2v_2)l (27) 
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FIG. 1: Schematic representation of the process T2. 



where n = (nL+n]i)/2 and the 'transmission coeiEcients' 
are given by 




((*«» 



FIG. 2: Current-noise ratio in units of V/T^ 
current as a function of applied magnetic field 
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for the spin 



Tla(^) = 



T2{lo) 



2 

1 



yrTTvrxAo 
[r+rx(r_r+ + A2)+r_ric^2]^ 

[T\T\ + {V_-T+fuj^]. (28) 



One can see that two kinds of processes contribute to the 
statistics. The coefficients Tis and Tia represent sym- 
metric and anti-symmetric tunneling events in which one 
unit of spin is transferred. Therefore, Tia controls the 
odd cumulants, in particular the average spin current 
whereas Tis only influences even cumulants like the noise. 
Odd cumulants obviously vanish for a symmetric system 
(r_ = 0) or zero magnetic field (A = 0). On the other 
hand, T2 involves voltage-independent processes which 
transfer two units of spin. Those are coupled spin-flip 
processes which do not involve electron tunneling and 
can only occur at finite temperatures, see Fig. [TJ 



V. SPECIAL REGIMES 
A. Zero temperature case 

Next, we shall address some special regimes and begin 
with the zero temperature case. In the limit of low energy 
we obtain the following statistics. 



1 



x(/i)= (1 - Ti,) + -Ti, ( 



,v 



-ili\ 



N 



(29) 



This is basically a symmetric trinomial distribution. 
Restoring SI units, the exponent can be written as 
N — eVT/h, in which case N can be interpreted as the 
number of 'attempts', ^^ i.e. the number of electrons im- 
pinging on the tunneling contact in time T. Comparing 
it with the corresponding distribution function for the 
charge transfer,— one concludes that ((29)) is the same dis- 
tribution with odd moments removed. Hence, the lowest 
order average spin current is quadratic in V, precisely as 
found by Schiller and Hershfield^' The noise in the zero 



temperature limit, however, is linear in V, 
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TV r+r_L (r_r. 



TT 


[A^ 


+ r+(r^ 
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1 - To] , 







(30) 



rT 



where Ss = /„ dtls{t) in the spin transferred during the 
measurement time T. Tq is an effective transmission co- 
efficient defined by To := ^^1^(0) + T2(0) which also 
emerges in the finite temperature case as we shall see 
soon. Equation pO|) is a manifestation of the general 
formula for the shot noise power of a tunneling junction 
with transmission coefficient Tq. The same leading or- 
der dependence on the applied voltage has already been 
obtained perturbativelyJ^ Moreover, in this limit, the 
same formula has already been obtained for the charge 
current fii 

Due to the different behavior of current and noise as 
a function of the applied voltage, one obtains a voltage- 
dependent current-noise ratio 



{{5s 



i/yrTTVA 



(31) 



{{Ss^)) r+(r_r+ + A2)- 

One notices that for constant magnetic field the ratio 
becomes more favorable for large bias voltages. The de- 
pendence on the magnetic field for constant bias voltage 
is plotted in Fig. [21 Notably, the function has a minimum 
at A = yTTr^. Towards large magnetic fields, the ratio 
vanishes. 

We can go one step further and calculate the third cu- 
mulant in the same limit, which is quadratic in V. Nor- 
malized with respect to the average spin current, 



{{Ss')) 



6- 



r+r^ ^_^^ 



A2 



1 



(32) 



{{Ss)) '[A^+T+{r^+T^)] 

which can be interpreted as a generalized Fano factor. 

B. Equilibrium case 

The other easily accessible regime is the equilibrium 
case V = 0. Here, we trivially expect all odd cumu- 



lants to vanish. Indeed, setting n^ = nji in (j27p . we 
notice that the anti-symmetric part Tia vanishes and the 
CGF becomes completely symmetric. The lowest order 
cumulant is therefore the second one for which one ob- 
tains exactly the same Johnson-Nyquist result as for the 
charge current fluctuations 



{{6s^)) = AGoTkeTTo, 



(33) 



where Go = e^ /h is the conductance quantum. By direct 
calculation one verifies that in spite of finite magnetic 
field the particle currents carrying different spins arc in- 
deed completely uncorrelated. 



VI. CROSS-CORRELATIONS 

Finally, we shall make use of the availability of the 
combined PCS to obtain the correlation between charge 
current and spin current. The quantity {{Ssdq}) can easily 
be obtained from ([M)) by double partial derivation with 
respect to /i and A. Again, for zero temperature and low 
energies, we obtain the normalized correlation 



{{SsSq)) 



^/WfW^lJ [A2 + r+(r_ + r^)f 



rir^ 



r_r. 



A^)' 



r+rx(r_r+ + A2) 



(34) 



This correlation function is plotted in Fig. [3l It can be- 
come negative and, in particular, changes its sign when- 
ever the sign of y or A is switched. This was expected as 
the substitution A ^ — A inverts the spin current while 
leaving the charge current unchanged whereas substitut- 
ing V ^ —V acts analogously on the charge current. 
Moreover, we notice that there is a choice of parameters, 
r+r^ = A2+r_r+, where this correlation becomes non- 
trivially zero. The correlation of spin current and charge 
current comes about as a result of the coupling of both 
currents to the dot level. Setting A = /i = in (fT4|) . we 
notice that such a coupling can be removed by setting 
r_ = A = 0. In this case, it can indeed be shown that 

X(A*,A) =Xs{f^)XqW- 

At T = 0, in linear response the combined statistics 
takes an appealing form. 



X(A, m) = [1 + Ti(cos^e''^ - 1) + T2ie^'^ - 1)] 



N 



where we recover the same transmission coefficients that 
had already been found in the case of charge transportii 
and which are responsible for electron and electron pair 
tunneling. 



Ti 

T2 



2r+ri(r_r4 
r^ri/D(o). 



A2)/D(0), 



(35) 



We can recover the familiar binomial distribution if we go 
over to a CGF for transfer of electrons with preselected 



({SsSq}} 
0.4- 




FIG. 3: Normalized spin-charge current correlation in units 
of V/r+ as a function of magnetic field for different values of 

r_ 



spin orientation. To that end we introduce a set of differ- 
ent counting fields with self-evident notation A|^x = A±^. 
Then in analogy to the reduction presented in 6lllll| we 
obtain 



x(At,x) - { [1 + T,(e^^T _ 1)] [1 + re(e'^i - 1)] } 



N 



where the effective transmission coefficient is defined as 
Te = \fTi- This trivialization, however, takes place in 
the linear response regime only. Nevertheless, this is a 
remarkable result since even in the spin-resolved case the 
statistics of electron transfer turns out to be determined 
by a single parameter Tg. 



VII. STRONG-COUPLING FIXED POINT 

At temperatures below the Kondo temperature Tr-, 
the dot spin hybridizes with the lead spins and forms 
a singlet state. At this so-called strong-coupling fixed 
point, the system can be described as a pure resonant 
level system. The free Hamiltonian can be written down 
in terms of non-interacting s- and p-wave electrons, c^a 
and ana respectivelyj-iS, 



^pa 



^0 / ^ ^pa \ppa^pf7 "^ ^pa^po 



V 



o / , [pprj^pc ' ^pa'^pc] 



(36) 



pa 



where ep^ ~ ep + ah includes the magnetic field h and V 
is the voltage difference between the two leads. The low- 
est order excitations away from the strong-coupling fixed 
point are given by a scattering term and an interaction 



term, 



Hf,r — 



Hi, 



2itvTk , 



/ A^P ~^ ^p/'^pG^P'' 



(37) 



pp' a 



tti/'^Tk 



(c|c| -no)(c|cx -no), (38) 



where no = (cJ.Cct)„ stands for the particle density at 
V = and h = 0. Moreover, i/ denotes the density 
of states and a and cf) are dimensionless aniphtudes. In 
the actual Kondo model the Fermi-liquid relation a — 
4> = 1 holds. ^^ Standard perturbation theory in a and (j) 
shows that no spin current is to be expected. Therefore, 
we will focus on higher moments which are conveniently 
accessible by means of full counting statistics. 

For this purpose we include the counting fields for 
charge and spin current. The most general approach in- 
volves the insertion of a spin-resolved charge counting 
field Act . Then, the approach of [18|[l9| can easily be gen- 
eralized to 

Co-(Acr) == cos(Ao-/4)ccr - isin(Acr/4)acr , 

ao-(Acr) = -isin(Ao-/4)ccr +cos(Ao-/4)acr . (39) 

The choice of this rotation is motivated by the fact that 
re-inserting the first order expansion in Ao- into the Hsc 
and Hint exactly gives the charge and spin currents pro- 
duced by these terms. We shall denote by 



Hx = H,,[c{\)]+H,nt[c[\)] 



(40) 



the perturbation Hamiltonian with rotated Ccr and a^ 
operators. One can easily see that this rotation leaves 
Ho invariant. 

The term Hq does not couple the spins and therefore, 
the unperturbed statistics is a trivial combination of the 
statistics for a spin-polarized system. Indeed, it is easy 
to show that 



Xo(A.)=e^^(^T+A.), 



(41) 



where N = (TF)/(27r) — eVT/h is the number of elec- 
trons of a certain spin impinging on the impurity during 
the measurement time T. In order to address the lowest- 
order corrections in a and (j), we start from the general 
expression 



x(A<t) =Xo(A<t) (Tcexp 



-i I dtHx{t) 
'c 



(42) 



A linked cluster expansion to the lowest contributing or- 
der in H\ leads to 

lnx(A,) = iN{\^+\{) 
1 



dtidt2{TcHx{h)Hx{t2)) . (43) 



c 



This expression contains terms proportional to a^, (jP' 
and a(j). As the scattering Hamiltonian Hgc does not 
couple up- and down-spins, the spin-resolved statistics is 
again a trivial combination of the two orientations. 



(5a In X = 



a^rV V^ + 4(7rT)2 + UK^ 
487rr|. sinh(y/2r) 



^ [(e-*^" - l)e^/2^ + (e^^- - l)e-^/2^' 



(44) 



By changing the basis of the counting fields and introducing counting fields for spin and charge current, ^ — (A|— A|)/2 
and A = (A| -I- A|)/2, respectively, one can readily perform the sum over a and one obtains 



(5q In X = 



247rT|. sinh(y/2T) 



{[cos(M)e-^^ - l]e^/2T _^ [cos(^)e^^ - l]e-^/2^} 



(45) 



Being even in /i, this expression will not produce any spin current but it will contribute to the spin current noise, 
albeit in a trivial way. More interesting effects can be expected from the interaction terms. Investigating the terms 
proportional to cj)'^, one finds that the whole contribution S^lnx = Inxi + l^Xa consists of an exchange term Inxi 
and a mean- field- like term In X2 ■ The latter is given by 



lnx2 = 



^Th^ 



V 



2t:TI sinh(y/2r) 



[[cosi^l)e-'^ - l]e^/2^ + [cos(^)e^^ - l]e-^/^^} 



(46) 



The mean-field-like contribution is mere combination of spin-polarized terms. However, the exchange term contains 
new pieces. 



Inxi = 



(j>^TV V^ + 4(7rT)2 
247rT|. sinh(y/2T) 

127rT|. sinh(y/T) 



{[cos(/i)e-*^ - l]e^/2^ + [cos(A*)e^^ - l]e-^/2^} 



(I 



'2i\ 



l]e 



V/T 



2iX 



l]e-^/^} 



e 



[cos(2^) - 1] . 



(47) 
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The first contribution describes single electron tunneling and has the same structure as the terms resulting from the 
scattering Hamiltonian. The second part can be interpreted as electron-pair tunneling. The fact that this term does 
not contain the spin counting field is not surprising. It accomodates doubled counting field which according to [a| can 
be ascribed to electron pair tunneling. It is natural to expect that the spins of this pair are oriented opposite to each 
other, so that the spin current fluctuations are suppressed. The last term of (j47l) is a thermal, voltage-independent 
contribution describing the transfer of two spins. Finally, the cross-term, which is proportional to acf), is given by 

S^^lnx = "!If ...„^,J.;.^^ {[cos(/i)e-^ - Ije^/^^ + [cos(M)e^^ - Ije-^/^^} . (48) 



.Ti 



sinh(y/2T) 



r 



VIII. COMPARISON OF THE RESULTS 

So far, we have analyzed the spin transport statistics 
for two distinct realizations of the Kondo system: (i) 
the Toulouse point corresponds to a certain quite special 
choice of the coupling constants Ji and (ii) the strong- 
coupling regime effectively describing the system in the 
case J — *■ oo, which, according to the RG analysis is 
reached at low energies. The purpose of this Section is 
to establish a connection between these two situations. 

There are several clear differences between the two 
regimes which lead to different transport statistics. First, 
the strong-coupling fixed point still possesses two sym- 
metries which are broken in the Toulouse fixed point by 
setting Jz^ ^ J^^ = and r_ ^ 0. Moreover, in 
the above calculations at the Toulouse point, the mag- 
netic field was only applied to the dot. This scheme can- 
not be easily implemented in the unitary limit since the 
dot level, being perfectly screened, effectively disappears 
from the system. Thus the magnetic field dependence in 
Section IVIII is referring to lead electrons. For this rea- 
son we shall only compare the statistics at zero magnetic 
field. 

From the Toulouse point calculation for A = r_ — 
0, we obtain the following expression for the combined 
statistics up to order V^, 



lnx(Ai,A) = 2^iVA+— -^ f 
dttI ^ ^ 



-2iA 



1 



(49) 



The corresponding strong-coupling result for h = Q reads 



lnx(/i,A) = 2iN\ + 
+ 



ja^ + 0^)Ty^ 
127rT|. 



[cos{fi) 



-iX 



6^r| 



-2iA 



-1] 

(50) 



Naturally, in the linear order in V these two terms are 
equivalent. Differences arise in the next-to-leading terms, 
most notable is the /x-dependence of (l50|) . The physical 
reason for this effect is the following: by discarding the 
J^^ and J^^ coupling terms at the Toulouse point, spin- 
flip tunneling is the only transport process available in 



the Toulouse limit. Under these conditions it is ener- 
getically favorable to tunnel electron pairs with opposite 
spins since this process effectively leaves the Kondo sin- 
glet state untouched. On the contrary, in the true unitary 
limit where Jz are finite, electron tunneling without 

the impurity spin-fiip is allowed, thus leading to a contri- 
bution containing a single A field which is also explicitly 
/i-dependent. Another insight which is worth mention- 
ing is the Kondo temperature estimate Tk ^ F^, which 
coincides with the original proposal of [8| . 



IX. CONCLUSION 

To summarize, we have investigated the combined 
statistics of charge and spin currents through a Kondo 
dot at the Toulouse point as well as in the unitary limit. 
In the Toulouse limit at zero temperature, the spin cur- 
rent noise grows linearly with V while the spin current 
is subleading. The spin transfer statistics is trinomial in 
the linear response regime. At finite temperatures, we 
have established the Johnson-Nyquist formula with an 
effective transmission coefficient for this type of trans- 
port. We have calculated explicitly the combined charge 
and spin current correlations, which turn out to change 
sign depending on system parameters. Furthermore, the 
spin-resolved charge transfer statistics is shown to be bi- 
nomial, independent on the spin orientation and governed 
by a single effective transmission coefficient. 

In the strong-coupling limit the spin-resolved FCS has 
helped to reveal a very rich structure of the transport pro- 
cesses. In addition to the expected spin-isotropic single 
electron tunneling, the contribution of electron pairs tun- 
neling in singlet spin configuration can be clearly made 
out. Moreover, the FCS contains a part entirely due to 
spin currents which surprisingly does not depend on the 
applied bias voltage and which vanishes at zero temper- 
ature. 
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